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ABSTRACT 

In this Letter we investigate uniformly rotating, homogeneous and axisym- 
metric relativistic fluid bodies with a toroidal shape. The corresponding field 
equations are solved by means of a multi-domain spectral method, which yields 
highly accurate numerical solutions. For a prescribed, sufficiently large ratio of 
inner to outer coordinate radius, the toroids exhibit a continuous transition to 
the extreme Kerr black hole. Otherwise, the most relativistic configuration ro¬ 
tates at the mass-shedding limit. For a given mass-density, there seems to be 
no bound to the gravitational mass as one approaches the black-hole limit and a 
radius ratio of unity. 

Subject headings: gravitation — relativity — stars: neutron — stars: rotation 
— black hole physics — methods: numerical 

1. Introduction 

Self-gravitating toroidal fluid configurations (without a central body) in Newtonian 
Gravity were the subject of analytic investigations by Poincare (1885a,b,c); Dyson (1892, 
1893); Kowalewsky (1895) and Lichtenstein (1933). In particular, Dyson was able to give a 
fourth order expansion of uniformly rotating, homogeneous and axisymmetric rings, which 
turned out to be an extremely good approximation for thin rings, see Ansorg et al. (2002b). 
Numerical evidence for the existence of these rings was given by Wong (1974) and later by 
Eriguchi & Sugimoto (1981), who, in addition, confirmed a conjecture by Bardeen (1971) 
stating that there is a continuous connection between the Maclaurin spheroids and the se¬ 
quence of ‘Dyson rings’, see also Ansorg et al. (2002b). 
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In this Letter we extend the above Dyson rings to Einsteinian Gravity. As one moves 
away from the Newtonian configurations, one observes typical relativistic effects such as the 
formation of ergo-regions. In particular, we found an interesting continuous transition to an 
extreme Kerr black hole. If a fixed ratio p\jpi > 0.5613 of inner to outer coordinate radius 
is prescribed and one gradually increases the gravitational mass (for fixed mass-density), the 
configurations eventually form an extreme Kerr black hole, as described by Mcinel (2002). 
If on the other hand pi/p 2 < 0.5613 is fixed, the ultimate configuration, as one increases the 
gravitational mass, rotates at the mass-shedding limit. It is furthermore interesting to note 
that the maximum gravitational mass seems to become infinite as p\ /p 2 —> 1. 

For computing the ‘relativistic Dyson rings’ we extended our multi-domain spectral 
method (Ansorg et al. 2002a) to toroidal topology. Again we obtain an accuracy of up to 12 
digits for configurations sufficiently far away from limiting cases. 

Relativistic Dyson rings and their possible generalizations (rings with realistic equations 
of state, configurations with a central object) might be relevant in two different astrophysical 
situations: (i) they could form as a result of stellar core-collapse in the case of high angular 
momentum, cf. Thorne (1998); (ii) they could be present in central regions of galaxies. Of 
course, for a final evaluation of the astronomical relevance of relativistic toroids their stability 
must be analyzed. It is to be expected that they are stable with respect to axisymmetric 
perturbations, but they may be unstable to (non-axisymmetric) fragmentation. In this 
case, the toroids could nevertheless play an interesting role as intermediate configurations in 
various astrophysical collapse scenarios. 

In what follows, units are used in which the speed of light as well as Newton’s constant 
of gravitation are equal to 1. 


2. Metric Tensor, Field Equations and Boundary Conditions 

For an axisymmetric and stationary space-time describing the gravitational field of a 
uniformly rotating perfect fluid body the line element can be cast into the following form 4 : 

ds 2 = e 2a (dp 2 + d( 2 ) + W 2 e~ 2u (dp — lo dt) 2 — e 2u dt 2 . 

We define the corresponding Lewis-Papapetrou coordinates (p, £,</?, t) uniquely by the re¬ 
quirement that the metric coefficients and their first derivatives be continuous at the surface 


4 This form of the line element of an axisymmetric and stationary space-time also applies in more general 
situations, e.g. for differentially rotating perfect fluids or if the stress-energy tensor is that of a stationary, 
axisymmetric electromagnetic field (Carter 1969). 
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of the body. 

A particular consequence of the interior held equations for a perfect fluid body revolving 
with the uniform angular velocity hi is the boundary condition 

e 2l/ — W 2 (u — fl) 2 e _2t/ = const. = (1 + Z 0 )~ 2 , 

which holds along the surface of the fluid. The constant is related as shown to the relative 
redshift Zq, measured at infinity, of photons which are emitted from the body’s surface and 
do not carry angular momentum. 

Interior and exterior held equations together with the above boundary and transition 
conditions at the body’s surface, asymptotic behaviour at infinity and regularity conditions 
along the rotation axis (p = 0) form a complete set of equations to be solved, see for example 
Butterworth & Ipser (1976). 

As long as we do not consider the transition to the extreme Kerr black hole, the regularity 
condition at infinity is asymptotic flatness. If we however follow the ‘physical route’ to the 
extreme Kerr black hole studied by Meinel (2002), we learn that in this limiting process the 
fluid body shrinks until it coincides with the coordinate origin. The geometry of the exterior 
space-time assumes that of the extreme Kerr solution outside the horizon. On the other hand, 
a completely different space-time, which is not asymptotically flat, forms if we rescale the 
coordinates such that the fluid body retains its finite extension. The asymptotic behaviour of 
the corresponding gravitational potentials in this limit is given by the ‘extreme Kerr throat 
geometry’, see Bardeen & Horowitz (1999) and Meinel (2002). In order to determine the 
physical parameters of the rings in the black-hole limit precisely, it is necessary to calculate 
this inner solution with its non-flat asymptotic behaviour. 


3. The Multi-Domain Spectral Method 

In the case of spheroidal figures of equilibrium, a two-domain spectral method was 
used to yield highly accurate numerical solutions (Ansorg et al. 2002a). In this method we 
separately mapped the interior and exterior of the star onto a square. The field quantities as 
well as the unknown shape of the fluid body were written in terms of Chebyshev expansions, 
and the corresponding coefficients resulted from a high-dimensional nonlinear set of equations 
that incorporated both field equations and transition conditions. This system was solved 
by a Newton-Raphson method, and an initial guess for the solution was taken from the 
analytically known Newtonian Maclaurin spheroids. 

The idea of mapping several subregions separately onto squares can also be applied to 
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the case of toroidal figures. In a first step we map the interior of the ring onto a square 
(0 < s < 1, 0 < t < 1) by 

P 2 = p\ + (pI-p\)s 
c 2 = (1 - t)y B (s) 

where pi /2 stands for the inner and outer coordinate radius of the ring, respectively, and the 
non-negative function y B , which describes the unknown shape of the ring’s surface, satisfies 


yn(0) — 0 , 2 /b(1) — 0. 


(We assume reflectional symmetry with respect to the plane ( = 0.) Next we introduce 
toroidal coordinates (p, () in order to obtain a compact coordinate region for the entire 
space exterior to the ring: 

z = ip m cot [z/ 2] 


with 


z = p + i£ and z = p + [(. 


The value p m must be chosen such that p\ < p m < p 2 - Since in these coordinates the metric 
potentials are not analytic at z — 0, it is necessary to divide the corresponding compact 
coordinate region of values (p, Q) into further subregions, each one of which is again to be 
mapped onto a square. 


The solution is represented and determined in a completely analogous manner as de¬ 
scribed above for the spheroidal bodies. Note that the initial guess for the Newton-Raphson 
method now comes from the numerically known Dyson rings in Newtonian Gravity. 

Again, we obtain very accurate solutions 5 yielding up to 12 digits for configurations 
that are sufficiently far away from the mass-shedding limit and from the limits pi/p 2 —*■ 0, 
Pl/P2 -»■ I- 


4. Results 

For a particular (constant) energy density /i, the relativistic Dyson rings are charac¬ 
terized by two parameters, say the redshift Zq and the radius ratio p\jp 2 - The region in 
which these parameters may vary is depicted in Fig. 1. Here, vanishing Z 0 represents the 
Dyson rings in Newtonian Gravity. The fluid bodies with pi/p 2 = 0 are transition con¬ 
figurations from toroidal to spheroidal topology. Starting from the Newtonian body, these 


5 The accuracy can be tested in several independent ways, see Ansorg et al. (2002a). 
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configurations reach a mass-shedding limit at Zq « 0.26. If we now follow the mass-shedding 
curve, Zq increases and reaches infinity at P 1 /P 2 = 0.5613, which corresponds to a transition 
to the extreme Kerr black hole. For a rigidly rotating disk of dust, such a transition was 
conjectured by Bardeen & Wagoner (1969, 1971) and analytically proven by Neugebauer & 
Meinel (1993, 1995). Also for differentially rotating disks of dust a transition of this kind has 
been found (Ansorg & Meinel 2000; Ansorg 2001). In this limit the coordinate extension of 
the gravitational source (here the relativistic Dyson ring) shrinks until the object coincides 
with the coordinate origin. The metric tensor assumes the form of the extreme Kerr solution 
outside the horizon 6 with W = p, and the angular momentum J, the gravitational mass M 
and the angular velocity D approach the relation 

J = M 2 = (2D) -2 . (1) 

This extreme Kerr black-hole limit emerges for radius ratios pi/P 2 > 0.5613. If we move from 
here along this boundary curve towards pi/p -2 = 1, we note that the normalized gravitational 
mass M = Mp 1 / 2 , the normalized rest mass 7 M 0 = Mop 1 / 2 and the relative binding energy 
(M 0 — M)/M 0 increase, see Fig. 2. In the Newtonian limit (Z 0 —> 0) the normalized angular 
velocity D = D/p 1 / 2 tends to zero as pi/p 2 —> 1. Our numerical results support the plausible 
assumption that the same is valid for all values of Zq. In the black-hole limit ( Zq —► oo) 
this would imply M —» oo as p\jP 2 —► 1, see Eq. (1). Therefore we believe that there is no 
bound to the gravitational mass of relativistic Dyson rings of a given mass-density as one 
approaches the black-hole limit and a radius ratio of unity. 

The evolution of a typical relativistic Dyson ring sequence with fixed radius ratio pi/p 2 = 
0.7 from the Newtonian to the black-hole limit can be seen in Figs 3 and 4. Corresponding 
values of various physical quantities are given in Table 1. Finally, in Fig. 5 we provide 
selected examples of configurations at the mass-shedding and the black-hole limit. More 
details of our methods and results including the discussion of realistic equations of state, 
configurations with a central object (black hole or neutron star), stability, and astrophysical 
relevance will be published elsewhere. 

This work was supported by the Deutsche Forschungsgemeinschaft (DFG-project ME 
1820/1-3). 


6 In our coordinates, the horizon (and the ‘throat’ of the extreme Kerr metric) is given by p = ( = 0, see 
Meinel (2002). 

7 For the calculation of M 0 we assume that the total energy density g is equal to the rest-nrass density, 
cf. Bardeen (1971). 
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Table 1: Physical quantities for the configurations with the radius ratio p\ /p 2 = 0.7 dis¬ 
played in Fig. 3. Here H = Q,/p 1 ^ 2 , M = Adp 1 ! 2 , M 0 = M 0 p 1 ^ 2 , J = Jp, and .R c i rc = R c \ rc p 1 ^ 2 
are normalized values of the angular velocity fi, gravitational mass Ad, rest mass Mo, angular 
momentum J, and circumferential radius R c i rc = We~ v [p — p 2 ,C — 0]. Note that M, M 0 , J, 
and R circ tend to zero in the Newtonian limit Z 0 —> 0, whereas 0, for pi/p 2 = 0.7, approaches 
the value 0 N = 0.48109. 


^0 

0 

M 

M 0 

J 

-f^circ 

0.05 

4.9108e—01 

7.9661e—03 

8.0842e—03 

2.3168e—04 

2.8483e—01 

0.50 

5.5491e—01 

1.4018e—01 

1.5659e—01 

2.8624e—02 

7.4643e—01 

1.22 

6.0702e—01 

2.8798e—01 

3.5082e—01 

9.8186e—02 

9.7512e—01 

1.60 

6.2320e—01 

3.3896e—01 

4.2564e—01 

1.3035e—01 

1.0441e+00 

2.50 

6.4726e—01 

4.2094e—01 

5.5610e—01 

1.9083e—01 

1.1543e+00 

6.00 

6.7932e—01 

5.5338e—01 

7.9987e—01 

3.1294e—01 

1.3466e+00 

13.0 

6.9211e—01 

6.2871e—01 

9.6211e—01 

3.9769e—01 

1.4734e+00 

oo 

6.9980e—01 

7.1449e—01 

1.1742e+00 

5.1050e—01 

1.6427e+00 
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Fig. 1.— Parameter region of the relativistic Dyson rings in the p\j P 2 - Zq/{1 + Zq) - plane. 
Rings with fixed radius ratio pi/P 2 = 0.7 are depicted in Fig. 3, see also Fig. 4 and Table 1. 
The boundary configurations indicated by «’s can be found in Fig. 5. 
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Fig. 3.— Cross-sections of relativistic Dyson rings with fixed coordinate-radius ratio 
P 1 /P 2 = 0.7 and varying redshift Z 0 . The normalized ^-coordinate DC is plotted against 
the normalized p-coordinate Dp (with the axes scaled identically). In the Newtonian limit 
(Zo —*■ 0) as well as in the black-holc limit (Z 0 —> oo) the ring shrinks down to the normalized 
coordinate origin. Dashed lines represent the boundary of the toroidal ergo-region. In the 
limit Zq —> oo the ergo-region is that of the extreme Kerr black hole, cf. Meinel (2002). 





Fig. 4.— The dimensionless quantities M 2 /J, OM and (Mo — M)/Mq for the relativistic 
Dyson ring sequence with pi/p 2 = 0.7 from the Newtonian limit (Z 0 —> 0) to the black-hole 
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Fig. 5.— Cross-sections of relativistic Dyson rings at the mass-shedding limit (a) - (c) and 
at the black-hole limit (c) - (e), cf. Fig. 1. Here the normalized ^-coordinate (/p 2 is plotted 
against the normalized p-coordinate p/p 2 (with the axes scaled identically). The dashed 
cross-section in (b) again represents the ergo-region. Note that Vtp 2 —► 0 in the black-hole 
limit, see Fig. 3. Hence, in contrast to (b), the ergo-region would become infinitely large 
with the scaling of this figure. It is interesting to observe a more and more circular shape 
of the ring’s cross-section as one approaches p\j p 2 — 1, a property known from the Dyson 
rings in Newtonian Gravity. 



